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Modeling line-driven disk wind for broad absorption lines of quasars 
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Abstract 

The disk wind, which is powered by the radiation force due to spectral lines (line force), is studied 
for broad absorption line (BAL) quasars. We investigate the structure of the disk wind based on the 
non-hydrodynamic method and compare with wind properties inferred from X-ray observations of BAL 
quasars. In this paper, we apply the stellar wind theory to the initial condition (the mass outflow rate at 
the base of the wind) . We found the funnel-shaped winds with a half opening angle of ~ 50° for the case 
of e = 0.3 — 0.9 and Mbh = 1O 7_8 5 M0, where e is the Eddington ratio and Mbh is the black hole mass. 
Thus, the absorption features are observed for an observer of which a viewing angle is around 50°. A 
probability of BAL quasars is ~ 7 — 11%, which is roughly consistent the abundance ratio of BAL quasars, 
~ 10 — 15%. Here, the probability is estimated by the solid angle, that the absorbing features would be 
detected, divided by An. In contrast, if the Eddington ratio is smaller than 0.01 or if the black hole is 
very massive, Mbh ^ 10 9 M Q , the disk wind is not launched due to the less effective line force. Then, the 
quasars are identified as non-BAL quasars independently of the observer's viewing angle. 

Key words: accretion: accretion disks — galaxies: active — methods: numerical — quasars: absorp- 
tion lines — radiative transfer 



1. Introduction 

Accretion disks surrounding massive black holes are ori- 
gins of activity of active galactic nuclei (AGNs). Their 
continuum emission comes from the disks and absorption 
as well as emission lines are thought to be produced by the 
matter above and/or around the disks. Broad absorption 
lines (BALs), of which line widths exceed 2,000kms _1 , 
are observed in ^ 10 — 20% of quasars (Wcymann ct al. 
1991; Hamann et al. 1993; Allen et al. 2011). The BALs 
are caused by metals moderately ionized and blue-shifted 
with typical speeds of ~ 10,000kms _1 , up to ~ 0.2c 
or more (e.g., Jannuzi et al. 1996). Weymann et al. 
(1991) suggested that the difference between the BAL 
and non-BAL quasars is caused by observers' viewing an- 
gles, since the properties of emission lines and continua 
of the BAL and non-BAL quasars are remarkably simi- 
lar. Elvis (2000) proposed a phenomenological model of 
funnel-shaped disk wind, in which the matter is blown 
away toward the direction of the polar angle of ~ 60 — 66°, 
to explain the abundance ratio of BAL quasars (~ 10%) 
via the viewing angle. If outflows launched from the sur- 
face of the accretion disks are accelerated toward the ob- 
server, the blueshift of the BALs could be explained. On 
the other hand, if the wind does not interrupt the line of 
sight, then BAL features would not emerge in the spectra. 



Many theoretical models have been proposed to account 
for the origin of the outflows so far. One plausible sce- 
nario is that of 'magnetically driven winds' (Blandford & 
Payne 1982; Konigl & Kartje 1994; Everett & Murray 
2007). In this model, the matter is magnetically acceler- 
ated. However, this model needs an extra mechanism to 
explain that the metals are moderately ionized, since the 
gas irradiated by the strong X-ray around the nucleus is 
fully photoionized (the so-called overionization problem) . 

Another plausible force that accelerates the disk wind 
is the radiation force due to spectral lines (line force). 
This model can explain both acceleration and ionization 
states. In this model, the gas on the surface of the ac- 
cretion disk is accelerated by absorbing ultraviolet (UV) 
radiation from the disk through the bound-bound tran- 
sition. Since the bound-bound absorption does not ef- 
fectively occur for the overionized metals, the line force 
can accelerate the gas only in the lower-ionization state. 
Indeed, Stevens & Kallman (1990) showed that the line 
force works efficiently for the gas in the lower-ionization 
state and could be much more powerful than the radiation 
force due to electron scattering. Thus, though the lumi- 
nosity of the disks in most quasars does not exceed the 
Eddington luminosity, the line-driven disk wind can be 
launched. The discovery of the line locking in the quasars 
(Foltz et al. 1987) strongly supports the idea that the line 
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force plays an important role in driving the outflows. 

Proga et al. (f998, f999) performed two-dimensional 
radiation-hydrodynamic (RHD) simulations of winds from 
disks around white dwarfs. This RHD method has been 
improved and applied to the disk winds in AGNs by 
Proga et al. (2000) and Proga & Kallman (2004) (here- 
after PK04) , in which the radiation transport of X-ray and 
UV is taken into consideration. They showed the funnel- 
shaped disk winds of which the opening angle is 6 ~ 70° 
where is the polar angle measured from the rotation 
axis. The line-driven wind is launched at the distance 
from the black hole of ~ several 100i?s, where i?s is the 
Schwarzschild radius, and the ejected matter goes away 
mainly in the direction of 9 ~ 70°. Their simulations re- 
vealed that the line-driven wind is in fact formed for a set 
of the black hole mass (Mbh = 1O 8 M0) and the Eddington 
ratio (e = 0.5). PK04 suggested that winds can be pro- 
duced for Mbh > 10 7 M Q , but no winds appear for e = 0.1. 
This parameter dependence should be explored for a wider 
range of the parameters. Schurch et al. (2009) and Sim et 
al. (2010) calculated the spectra and compare them with 
those of the X-ray observations. However, they also only 
investigated the case of e = 0.5 and Mbh = 1O 8 M . 

Risaliti & Elvis (2010) studied the disk winds in AGNs 
for the wide parameter space of the black hole mass and 
the Eddington ratio. They investigated steady structures 
by non-hydrodynamic calculations, in which they solved 
the trajectories of the matter ejected from the disk sur- 
face without calculating the RHD equations. Their results 
seem to be consistent with the results of the RHD simula- 
tions. This non-hydrodynamic method is a powerful tool 
to investigate the dependence on several unknown param- 
eters. Here, we use a similar method to investigate BAL 
quasars. However, Risaliti & Elvis (2010) did not research 
the wind properties, the ionization parameter, the veloc- 
ity and the column density along the line of sight. Thus, 
within the framework of the non-hydrodynamic model, it 
is worth studying whether the disk wind explains the BAL 
features. 

In Risaliti & Elvis (2010) and PK04, the density and 
the velocity at the base of the wind (at the disk surface) 
are treated as free parameters and simply assumed to be 
independent of the distance from the black hole. However, 
the stellar wind driven by the line force was investigated 
in detail. CAK75 derived the mass outflow rate (the den- 
sity and the velocity) as functions of the gravity and the 
radiative flux. Thus, we should employ the mass outflow 
rate of the wind base consistent with the prediction of 
CAK75. Based on the CAK theory, the density at the 
wind base is not constant but decrease with distance from 
the black hole, if the velocity is assumed to be the sound 
velocity (see section 2.2). 

In the present paper, employing the wide range of the 
parameters of the black hole mass, the Eddington ratio, 
we investigate the conditions under which the disk wind 
model can reproduce the X-ray absorption features of the 
BAL quasars. Although we do not study the spectra, we 
calculate the ionization parameter, the velocity and the 
column density for wide e-Mbh range and compare them 



with those inferred by the X-ray observation. 

We describe our calculation method in § 2. We present 
results in § 3. Finally, § 4 and § 5 are devoted to discussion 
and conclusions. 

2. Numerical Methods and Models 

2.1. Outline 

We study a steady structure of line-driven disk winds 
and apply them to the BAL quasars. We investigate the 
conditions under which the disk wind model can explain 
the X-ray absorption features seen in the BAL quasars. 

The method for investigating the structure of the wind 
is basically the same used by Risaliti & Elvis (2010), but, 
based on the stellar wind theory, we employ a more re- 
alistic mass outflow rate proposed by Castor, Abbott, & 
Klein (1975) (hereafter CAK75). Furthermore, we per- 
form higher resolution calculations and carefully treat the 
gas temperature and the ionization state to compare the 
model with wind properties inferred from the observa- 
tions. 

The steady structure of the disk winds is given by the 
calculation of the trajectories of fluid elements (stream- 
lines), which are ejected from the surface of the geomet- 
rically thin and optically thick disks (Shakura & Sunyaev 
1973). We solve the equation of motion considering the 
line force, coupling with the mass conservation along the 
streamline. We also calculate the temperature and the 
ionization parameter by taking account of the radiation 
transport of the X-ray from the vicinity of the black hole 
and UV emitted at the inner part of the disk. We show 
the schematic picture for our method in Figure 1. We 
calculate the streamlines in order of the distance from the 
black hole. When we solve i-th streamline, we take into 
consideration the extinction of the radiation via every in- 
ner streamline of j-th (j < i) as well as the self-shielding by 
i-th flow. We can obtain the global structure of the disk 
winds after calculating the final (outermost) streamline. 

Next, we calculate the column density, the outward 
velocity, and the ionization parameter for a wide range 
of viewing angles and compare them with those inferred 
from X-ray observation of the BAL quasars. We evaluate 
the range of the viewing angle (i.e., the solid angle), at 
which our wind model is consistent with the observations. 
The solid angle divided by 47r means the probability that 
the system is identified as a BAL quasar (hereafter, BAL 
probability). We investigate the BAL probability for a 
variety of the black-hole mass and the Eddington ratio. 

2.2. Initial conditions and basic equations 

By solving the equation of motion including the radi- 
ation force, we investigate the trajectories (streamlines) 
of the fluid elements ejected from the disk surface in the 
cylindrical coordinates, (r,ip,z). Here, the black holes and 
the accretion disk are located at the origin and the equa- 
torial plane (z = 0). We assume the system to be axisym- 
mctric with respect to the rotation axis. 

The initial position (wind-base) of a fluid element is 
(^O)^o)- The first streamline starts at the innermost ra- 
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Fig. 1. Schematic picture for our calculation method. The 
black hole and the accretion disk are located at the origin and 
the equatorial plane (z = 0). Solid orange lines indicate the 
streamlines, and AL is the geometrical thickness of the each 
flow. The length of the light ray passing through the j'-th flow 
is Asj. When we calculate the i-th streamline, the extinction 
of the radiation by inner streamlines (1st, 2nd,... i — 1th) and by 
the self-shielding of the i-th flow is taken inoto consideration. 



dius of tq = 20i?s (Here, we note that our results do not 
change so much if we employ ro = 10i?s). We calculate 81 
streamlines by setting the interval of the initial points in 
r-direction to be Ar = 8i?s (We confirmed that the wind 
structure does not change significantly if we set the inter- 
val to be Aro = 2i?g). That is, the initial radius of the z-th 
streamline is 20i?s + 8(i — l)Rs, an d the final fluid element 
is launched at the radius of ro = 660i?s- The initial alti- 
tude, zq, is set based on the hydrostatic balance in the ver- 
tical direction between the gravity and the pressure forces 
of gas and radiation, GMbuZq/tq = c s /zq + (a c /c)aT^ s , 
where G is the gravitational constant, c s is the sound 
speed of the gas at the mid-plane (z = 0) , a e is the mass- 
scattering coefficient for free electrons, c is the speed of 
light, (7 is the Stcfan-Boltzmann coefficient, and T e g is 
the effective temperature. In the top panel of Figure 2, 
we show zq as a function of the distance from the black 
hole, r, normalized by i?s being the Schwarzschild radius, 
where we employ the alpha viscosity parameter of 0.1. We 
find zq depends on the black hole mass and the Eddington 
ratio (mass accretion rate). 

At the wind base (r = ro, z = zq), we set the horizon- 
tal and azimuthal components of the velocity to be null 
and the Kcplerian velocity. The initial density, po, and 
the vertical component of the initial velocity, vq, are set 
so as to meet the prescription by CAK75, in which they 
investigated the steady and spherical wind driven by line 
force. The resulting mass outflow rate per unit surface at 
the distance of R is 

1 GM* \(l-cO/a 



a c v th R 2 

x(fcr) 1 / a (i 



r) 



-(l-a)/a 



(1) 



where Vth is the thermal speed of the gas, M„ is the mass of 
the star, a(~ 0.6) and 0.03) are the constants related 
to the line force, V is the Eddington ratio of the star [see 
equation (21) of CAK75]. By making minor revisions, we 
apply above relation to the mass outflow rate from the disk 
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Fig. 2. The radial profile of the initial altitude (top panel), 
the initial density (middle panel), and the initial upward ve- 
locity (bottom panel) are ploted for e = 0.3 and Mbh = 10 7 Mq 
(red line), for e = 0.5 and M B h = 1O 7 M (blue line), and for 
e = 0.3 and M B h = 10 s M (green line). 

surface. The gravity GM^/R 2 is replaced by GM Bii z /r 3 , 
which is the gravity at the wind base. We use the local 
Eddignton ratio, 



r' = 



cr e O-T e 4 ff /c 



3e 



GMbhzq/ R 3 4 V z /R s 7 
as substitute for E As a result, we have 
1 GMbhZq 



(2) 



Povo 



x(fcr') 1/Q (l 



a{l-a) {l - a)/a 

p' \-(l-a)/a 



(3) 



Throurghout the present study, we assume the vertical 
component of the initial velocity to be the sound veloc- 
ity at the disk surface, vo = (/cH^eff / /wip) 1 1 2 , where fcs is 
the Boltzmann constant, /i(= 0.5) is the mean molecular 
weight, and m p is the proton mass. Thus, the initial den- 
sity, po, is given by the equation (3). In Figure 2, we show 
Po (middle panel) and v (bottom panel) as a function of 
r/Rs- Since vq is much smaller than the escape velocity, 
the disk wind appears if the fluid elements are accelerated 
by the line force. 

Here we note that the density and velocity of the wind 
base are still unknown, although we apply the CAK75 
relation. For instance, the mass outflow rate might be 
enhanced if the magnetic pressure force cooperates for 
launching the disk wind. High resolution radiation (mag- 
neto) hydrodynamic simulations would make clear the 
point. 

The equation of motion including the radiation force is 



dv 



GM BH R 



,7T 



cr c (F uv + Fx) 



M fe£uv (4) 



R 3 r" 1 c c 

where v — (v T ,V v ,v z ) is the velocity, R = [(r 2 + z 2 ) 1 / 2 ] is 
the distance from the origin, R is the point vector of the 
fluid element, I is the specific angular momentum, Fx and 
Ftjv are the X-ray flux and the UV flux, and M is the 
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force multiplier (see below). The radiation force due to 
the free electrons is the third term of the right-hand side. 
The final term is line-driven radiation force, which may 
cause the disk wind in our model. The gas pressure is 
neglected. 

Here, we speculate that the wind flow is launched from 
the disk surface perpendicularly, and is then directed away 
from the central source. When the flow is perpendicular 
to the disk plane, the geometrical thickness of each flow 
(AL, see Figure 1), which is measured at perpendicular 
to the streamline on r-z plane, is kept constant, AL = 
Aro ■ If we simply assume AL oc r when the wind is blown 
outward, the relation of AL = (r/ro)Aro can describe the 
geometrical thickness of the flow both in the vertical and 
outward phases. Then, the density is given from the mass 
conservation equation (2nrALp\v\ = const) as 



is calculated by 



p(v 2 r + v 2 z ) 1/2 r 2 = const. 



(5) 



For comparison, wc employ the other model, in which 
we assume the width of the flow docs not broaden (the 
width is constant) not only when the flow is vertical 
to the disk, but also when the flow is blown outward; 
AL = Aro- ln this case, the density is estimated by 
p(v 2 + v 2 ) 1 / 2 r = const. In the latter model, the density 
of the wind matter is relatively large at the distant region 
(r » /'o). Thus, the maximum of the ionization param- 
eter becomes smaller. However, since the large fraction 
of the moderately ionized region appears where the flow 
is vertical to the disk (r ~ ro), the BAL probability (see 
the §2.5) is nearly independent of the model of the wind 
geometry 

2.3. Radiative fluxes of X-ray and Ultraviolet 

We consider two spectral components of the radia- 
tion, X-ray and UV. The luminosities of X-ray (Lx) 
and UV (i/uv) are given by Lx = /xeiedd and £uv = 
(l-/x)e£edd- Here, / x is the X-ratio, L x /(Lx + Lu V ), e is 
defined as e = (Lx + L\jy) / L e dd, and £ e dd is the Eddington 
luminosity. In the present study, e is a parameter and we 
employ fx = 0.15. 

X-ray is assumed to be emitted near the central black 
hole and treated as a central point source. The radial 
component (i?-component) of the X-ray flux is written as 

^ = ^e- (6) 

where tx is the optical depth for X-ray measured from 
the origin. The other components of the X-ray flux are 
zero. We consider that the radiation is attenuated by the 
wind flow in the present study. When we calculate the 
i-th streamline, we take into consideration the obscura- 
tion by every inner streamline of j-th (j < i). That is, 
if a ray from the center to the fluid element gets across 
the j-th streamline (j < i), the radiation suffers from the 
dilution. Even if the ray does not get across, the j-th flow 
sometimes contributes to the obscuration because the flow 
has geometrical thickness (see § 2.2). We also consider the 
self-shielding effect, which is extinction caused by the up- 
per stream of the i-th flow. Therefore, the optical depth 



E Arx - 



(7) 



j<i 



where Arx.j = paxAsj , with ax being the mass extinction 
coefficient for X-ray and Asj being the length that the ray 
passing through the j-th flow (see Figure 1). The mass 
extinction coefficient is defined by ax = a c for £ > 10 5 and 
ax = lOOtTo for £ < 10 5 , where £ is the ionization parameter 
defined as, 



£ [erg cms 



pR 2 



(8) 



This simple treatment of ax reflects that the photoelec- 
tronic absorption is not effective when the ionization pa- 
rameter is very high, since the number of bound electrons 
is too small (see also Risaliti & Elvis 2010). We have the 
r- and z-components of the X-ray flux is F x = F^(r/R) 
and F x = Fx( z l R)- The azimuthal component is zero. 

We suppose that an optically thick and geometrically 
thin accretion disk is the only UV source where the effec- 
tive temperature follows the radial profile of 



T e s = Tj„ ( — 



-3/4 



(9) 



where r- ln (= 3i?s) is the disk inner radius and T; n is the 
effective temperature at r = r; n . Wc set T; n so as to meet 
the condition of 



(l-/x)ei, 



dd 



2TrraT oS dr, 



(10) 



and the UV flux from the disk in the optically thin media 
is calculated by 



-Ft.Tii 



thin 



^ndn, 



(11) 



where a is the Stcfan-Boltzmann coefficient, n is the unit 
vector, and f2 is the solid angle. Here, note that wc inte- 
grate the radiation of the disk within r = r out , at which 
the effective temperature of the disk is 10 4 K, since the 
outer part is too cold to emit UV photons effectively, and 
since the disk luminosity of the outer part is much smaller 
than that of the inner part. For simplicity, we attenuate 
the radial component (i?-component) of the UV flux using 
the optical depth measured form the origin, t\jv, as 

(12) 



'thin 6 



The numerical method for calculating t\jv is the same as 
that for X-ray optical depth, but we use a e instead of ax- 
We assume that the dilution of the polar component of 
the UV flux is negligible, 

^UV = F thim (13) 

since the wind is expected to flow near the disk surface. 
The azimuthal component of the flux is zero by definition. 
In numerically, we divide the surface of the accretion disk 
into 256 x 256 small surface elements and calculate the 
UV flux. 
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2.4- Force multiplier 

In order to estimate the line force, we employ the force 
multiplier proposed by Stevens & Kallman (1990), which 
is modified version of CAK75. The force multiplier, which 
is ratio of line force to electron-scattering force, is written 
as 



M(t,0 = kt 



-0.6 



(1+^max) 



0.4 



1 



(tn I -' 

V^'/max / 



Here, t is the local optical depth parameter 
dv -i 



t = <j c pv th 



(14) 



(15) 



logi ?7n 



(16) 



,(17) 



where Vth = (ksT '/ fimp) 1 ^ 2 is thermal velocity with T be- 
ing the gas temperature, and dv/ds is the velocity gradient 
along the light-ray. Since k and ?y max depend on £ as, 

k = 0.03 + 0.385 exp(-1.4£ - 6 ), 

and 

6.9exp(0.16£ a4 ) log£<0.5 

9.1exp(-7.96 x 10~ 3 £) log£ > 0.5 

the force multiplier is the function of the local optical 
depth parameter, i, and the ionization parameter, £. 
In our calculation, we approximate the velocity gradient 

dv _dv x 

d^-dP (18) 
where dv\/dl is the velocity gradient along the stream- 
lines. As we will show in Figure 3, the fluid element is 
lifted up in nearly vertical direction, since the vertical ra- 
diation flux from around the wind base mainly accelerate 
the matter upwards. Thus, the velocity gradient along 
the right-ray would be approximated by that along the 
streamlines. Subsequently, the matter is blown away in 
the radial direction. Then, the inner part of the accre- 
tion disk mainly illuminate the outflowing matter. Thus, 
above approximation would be valid. 

We calculate the gas temperature assuming that the gas 
is in the radiative equilibrium, since the radiative cool- 
ing/heating timescale is much smaller than the cooling 
timescale via the adiabatic expansion (Proga et al. 2000), 

™ 2 (Gcompton + G X -£b,l) = 0, (19) 

where, Gcompton is Compton heating/cooling rate, 

G C ompton = 8.9xl0- 36 £(T x -4T), (20) 

Gx is the rate of X-ray photoionization heating and re- 
combination cooling, 

G x = 1.5 x 10~ 21 £ 1/4 T- 1/2 (1 - r/T x ), 



.Lb,i is the bremsstrahlung and line cooling, 
L h i = 3.3 x 10- 27 T X / 2 + 1.7 x iq-^-^t- 1/2 



xexp(-1.3 x 10 5 /T) + 10 



-24 



(21) 



(22) 



We assume the temperature of the X-ray radiation, Tx = 
10 8 K, which is the same as Proga et al. (2000). Since this 
method for calculating the gas temperature is almost valid 



in a low-density and high-ionization parameter regime, the 
gas temperature would be underestimated in the region of 
the high-density and low-ionization parameters. Thus, we 
set the lower limit of the gas temperature as being equal 
to the disk effective temperature, T(r,z) m - m = T c ff(r). 

Although, in the present work, the force multiplier de- 
pends on the gas temperature through the sound speed 
[see equation (14) and (15)], another model for the line 
force, which does not depends on the gas temperature, 
has been proposed by Gayley (1995). Our results does 
not change so mach if we employ another force multiplier, 
since the gas temperature is kept around 10 4 K in the 
region that the matter is mainly accelerated. 

2.5. Assessment of BAL probability 

A probability for observing BAL (BAL probability) is 
evaluated as follows. Using the resulting structure of the 
disk wind, we investigate the ionization parameter, the 
outward velocity, and the column density along lines of 
sight (from an observer to the central black hole) . The X- 
ray observations of the BAL quasars reported the ioniza- 
tion parameter to be < 500 (Reeves et al. 2003, Gallagher 
et al. 2004) and several x 1000 (Braito et al. 2006, Wang 
et al. 2008). The X-ray observations also reveal that the 
outward velocity and the column density of the gas in such 
lower-ionization state are larger than 10 4 kms _1 and larger 
than 10 23 cm -2 (Chartas 2007 and references therein). 
Thus, we consider two conditions: (A)the outward ve- 
locity of the matter with £ < 100 exceeds 10 4 kms _1 and 
(B)the column density of the gas with £ < 100 is larger 
than 10 23 cm -2 . We also add the condition that (C)the 
column density is smaller than 1.5 x 10 24 cm -2 , since the 
Compton thick objects would not be identified as BAL 
quasars. If three conditions arc satisfied, we suppose that 
the X-ray absorption features emerge in the spectra. If 
this is the case, our wind model could also explain the 
BALs in UV band, since the X-ray absorption closely cor- 
relate to the BALs seen in UV wavelength (Brandt et al. 
2000). We evaluate the solid angle, ^bal, in which above 
conditions are satisfied. The BAL probability is given by 
^bal/4tt. 

3. Results 

3.1. Structure of line- driven wind 

In Figure 3, red and black solid lines show trajectories 
of fluid elements (streamlines) ejected from the surface of 
the accretion disk in r-z plane. The mass of the black hole 
and the Eddington ratio, which are parameters in our cal- 
culation, are given as M B h = 10 7 Af Q and e = 0.3 (hereafter, 
we call these 'baseline parameters'). The initial density, 
velocity, and altitude of streamlines are shown in Figure 
2 (red lines). On the red lines, the matter is accelerated 
and attain a velocity equal to the escape velocity, v = v csc 
where v = (u 2 + v 2 + v 2 ) 1 / 2 , at the points of the red filled 
circles. In contrast, the flow velocity does not reach the 
escape velocity on the black lines. 

In the inner region, r < 30i?s, it is found that the ejected 
fluid element immediately returns to the disk surface (we 
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Table 1. The BAL probability (%) 
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Fig. 3. The steady structure of the line-driven disk wind for 
e = 0.3 and Mbh = 10 7 Mq. The black hole and the accre- 
tion disk are located at the origin and the equatorial plane 
(z = 0). The red solid lines show trajectories of fluid elements 
(streamlines) ejected from the surface of the accretion disk. 
The red circles on the red lines indicate the point where the 
wind velocity exceeds the escape velocity. At the inner region 
°f r IS 30-Rs> the ejected gas immediately goes back to the 
disk surface (black lines). It is also found that the wind is not 
launched at the outer region, r > 350i?s i although we can not 
see the streamlines in this figure. The BAL are observed at 
the viewing angle of 9 = 45.6° — 54.0° (between dotted-lines) . 



call this 'failed wind') (see black lines). In this region, the 
initial density is too high for the line force to lift up the 
matter. High density reduces the ionization parameter 
because of £ oc p~ x [see equation (8)]. Indeed, we have £ is 
much less than 100 in this region. In the lower ionization 
regime (£ < 100), the force multiplier decreases with an 
increase of the density, M oc p~ 06 , since M oc t~°- 6 where 
t oc p [see the equation (15)]. Therefore, the line force 
does not exceed the gravity and the gas is not blown away. 
Note that the mechanism of the failed wind in our work 
is different from that in previous works (Murray et al. 
1995; Risaliti & Elvis 2010). We will discuss the point in 
section 4.2. 

In the middle region, 30i?s < r < 350i?s, we find that the 
disk wind is successfully launched. The matter launched 
from the disk surface is lifted up almost perpendicular to 
the disk plane, and exceeds the escape velocity (filled cir- 
cles). Subsequently, such matter is blown away towards 
the direction with polar angle of 9 ~ 45° — 55°. The wind 
structure is understood as follows. Since the initial den- 
sity (Po ~ 10~ 14 — 10 -13 ) is smaller in this region than in 
the inner region, and since the ionization parameter still 
less than 100, the force multiplier is large (M oc p -0 6 ), 
typically several xlOO. Thus, the vertical component of 
the radiation force (line force) is stronger than the that 
of the gravity. In contrast, the radiation from the vicinity 
of the black hole is attenuated by the failed wind, thus 
the radial component of the radiation force can not play 
an important role. As a consequence, the matter is accel- 
erated in the vertical direction. Such a situation changes 
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at the region of 9 > 55°. The radiation from the inner 
accretion disk, which works to accelerate the matter out- 
ward, increase with a decrease of the polar angle and is 
not attenuated by the failed wind in the region. Therefore, 
outward radiation force bend the streamlines to radial di- 
rection and producing the funnel-shaped wind. 

In the outer region beyond r ~ 350i?s, the temperature 
of the accretion disk is less than 10 4 K and the disk around 
the wind base does not effectively emit UV photons. In 
addition, the UV radiation from the inner region of the 
accretion disk is obscured by the failed wind (r < 30i?s) 
and the disk wind (30i? s < r < 350i? s )- Thus, although 
the initial density as well as the ionization parameter is 
small, the line force is made powerless and fails to form 
the wind. 

According to the method mentioned in § 2.5, we inves- 
tigate viewing angles for which the X-ray absorption fea- 
tures of the BAL quasars are observed. In the case of the 
disk wind with the baseline parameters, the conditions 
(A), (B), and (C) arc satisfied at the viewing angle from 
9 = 45.6° to 54.0° (see §3.1). That is, then, the absorp- 
tion features are observed, since the lower-ionization mat- 
ter with sufficient radial velocity obscures the nucleus (see 
Figure 3). The BAL probability, £^ral / ^ ■> i s obtained as 
cos(45.6°) - cos(54.0°) - 11%. 

For the viewing angle smaller than 45.6°, there is no 
wind between the observer and the nucleus. When the 
viewing angle is between 54.0° and 63.0°, the column 
density is too large to detect the absorption features 
(Ah > 1.5 x 10 24 cm~ 2 ), since the failed wind, of which the 
density is very high, obscures the nucleus. Thus, although 
the conditions (A) and (B) are satisfied, the object is not 
identified as the BAL quasars. For large viewing angle , 
9 > 63.0°, neither the condition (A) nor (C) is satisfied. 

3.2. BAL Probability 

Table 1 presents the BAL probability. For e = 0.3 - 0.9 
and Af B H = 10 7 ~ 8 - 5 M Q , we find that the BAL probabilities 
are 5 — 11% and are not so sensitive to the black hole mass 
as well as the Eddington ratio. On the other hand, the 
BAL quasars are not observed (BAL probability is null) 
in the case of M BH > 1O 9 M or e < 0.01. 
3.2.1. Eddington ratio dependence 

For e = 0.01, the radiation is too weak to accelerate the 
matter. Since the wind is not successfully launched, the 
absorption features are not observed. The reason why the 
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Fig. 4. Same as Figure 3, but for e = 0.5 (left panel). The disk structure is almost same with Figure 3, although the larger 
Eddington ratio is employed. The right panel is also the same as Figure 3, but we employ the initial altitude for e = 0.5 [case (a)], 
the initial density for e = 0.5 [case (b)], and the disk luminosity for e = 0.5 [case (c)]. The larger initial altitude and density work to 
increase the opning angle of the funnel-shaped wind. In contrast, the larger luminosity tends to accelerate the matter upward. 



BAL probability is roughly kept constant for e = 0.3 — 0.9 
and Mbh = 10 7 ~ 8 ' 5 M© is understood by the left panel in 
Figure 4. This panel is the same as the Figure 3, but for 
e = 0.5. It is found that the wind launching region broaden 
as an increase with s. However, we find that the opening 
angle and the thickness of the funnel-shaped wind does 
not change so much. Therefore, the BAL probability is 
roughly independent of the Eddington ratio. 

Since the initial values and the disk luminosity changes 
with e, here, we at first consider an influence of the ini- 
tial density on the wind structure. For this purpose, we 
employ the initial density profiles for e = 0.5 (blue line 
of the middle panel in Figure 2), but other initial con- 
ditions and the disk luminosity are set to the values for 
s = 0.3 (red lines of the top and bottom panels in Figure 
2). The result is plotted by purple lines in the right panel 
of Figure 4 [case (a)]. We find that the resulting disk wind 
has an wide opening angle {9 ~ 65°). Since the initial den- 
sity increases as an increase of e (see Figure 2), and since 
the line force is less effective for the gas with large density 
(Mocp -0 ' 6 , see §3.1), the height of the failed wind is small 
and the stream lines of the disk wind bends at the small 
altitude for the case (a) in comparison with the case with 
baseline parameter. We conclude that the higher density 
tends to make the opening angle of the disk wind larger. 

Next, we consider an effect of the initial altitude, Zq. 
The blue lines of the right panel of Figure 4 show the re- 
sult for the case (b) in which the initial altitude for e = 0.5 
(blue line of the top panel in Figure 2) is supposed, al- 
though the luminosity, initial density and velocity remain 
the values for e = 0.3 (red lines of middle and bottom pan- 
els in Figure 2). That is, the larger zq is employed in com- 
parison with the case with the baseline parameters. Then, 
since the initial position is distant from the disk surface, 
the radiation force becomes less effective relatively and 
the disk wind with the large opening angle is generated. 
As we have discussed above, the less effective radiation 



force works to form the lower failed wind, leading to the 
disk wind with larger opening angle. 

The large disk luminosity tends to narrow the open- 
ing angle. The green lines [case (c)] show the stream- 
lines when the radiation fields are calculated by assuming 
e = 0.5 and other initial conditions remain the values for 
e = 0.3 (red lines in Figure 2). Comparing with the stream- 
lines for e = 0.3 (Figure 3), the matter, which is launch- 
ing at r < 150i?Sj is blown away more upward, inducing 
the small opening angle of the funnel-shaped wind. This 
is because the UV radiation is strong and the line force 
effectively accelerates the matter. We also find that the 
wind is launching from wide region of the disk, r < 400i?s, 
since the UV emitting region (T e g > 10 4 K) expands with 
an increase of e. 

To sum up, the increase of the luminosity, which con- 
tributes to reduce the opening angle of the wind, counter- 
acts the increase of the initial density and altitude, which 
works to expand the funnel. As a result, the BAL proba- 
bility is insensitive to the Eddington ratio for e = 0.3 — 0.9 
and Mbh = 1O 7_85 M0. Here, we note that the initial ve- 
locity has little influence on the wind structure. Indeed, 
we verified that the wind structure does not change if the 
initial velocity for e = 0.5 — 0.9 is employed. 
3.2.2. Black hole mass dependence 

The BAL probability is insensitive to the black hole 
mass for Mbh = 10 7 ~ 8 ' 5 M Q , but tends to be null for 
Mbh ^ 10 9 M Q . Figure 5 shows the streamlines for Mbh = 
10 s M Q and e = 0.5. This figure shows that the open- 
ing angle for Mbh = 10 8 M Q is comparable to or slightly 
smaller than that for Mbh = 10 7 Af Q . We also find that 
the launching region, which is measured in an unit of i?s , 
shrinks with an increase of the black hole mass. Such be- 
havior is caused by the decrease of the force multiplier as 
well as the temperature of the disk. 

The local optical depth parameter, t = o c Vthp{dvi / dl) -1 , 
is roughly evaluated as t — cr c poZo near the wind base, 
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Fig. 5. Same as the left panel of Figure 4, but for 
Mbh = 10 8 Mq . The half opening angle of the wind does not 
depend on the black hole mass. However, the wind launching 
region shrinks with an increase of the black hole mass. 

where we suppose Dth ~ w 0j P ~ Po- an d dvi/dl ~ vq/zq- 
Since zq is roughly proportional to Mbh {zq k Rs), T' is al- 
most independent form the black hole mass (see equation 
2). We have p oc (1/vq) (M| H /r 3 ) ocMg^ 4 with using the 
Thus, the optical depth parameter 



relation of vq oc M^/ S 



is proportional to and we have a very weak Mbh- 



-0.6 



oc M, 



-0.15 
BH 



dependence of the force multiplier, M oc f 
As we have discussed above, the disk wind with the wide 
opening angle form since the force multiplier is small (ra- 
diation force is less effective). Thus, the opening angle 
for Mbh = 1O 8 M is comparable to or slightly larger than 
that for Mbh = 10 7 M o . 

The size of the UV emitting region of the accretion disk 
(T c g > 10 4 K), which is normalized by i?s, is proportional 

— 1/3 

to M BH ' , in the case that the Eddington ratio is kept 
constant. This is the main reason that the launching re- 
gion of the wind becomes small when the black hole is 
massive. 

Since the UV emitting region becomes narrowed, and 
since the failed wind region slightly expands via the re- 
duction of the line force (M oc Mg^' 15 ), the disk wind 
tends to disappear when the black hole is too massive 
(Mbh ^ 10 9 M©). Hence, it is implied that the absorp- 
tion features are not identified in the quasars with very 
massive black holes. 

4. Discussion 

4-1. Comparison with observations 

Our model indicates that the BAL probability is several 
to 10 percent for a wide range of the Eddington ratio and 
the black hole mass, e = 0.3 - 0.9 and Mbh = 1O 7 " 8 ' 5 M . 
This result is roughly consistent with the abundance ratio 
of BAL quasars, - 10 - 20% (e.g., Allen et al. 2011). The 
resulting BAL probability might increase, if the nucleus 
is obscured by the dusty torus, which aligns with the ac- 



cretion disk, is not identified as quasars. For instance, 
the BAL probability is doubled, 10 — 20%, if the observer 
with 8 > 60° cannot detect quasars via the obscuration 
by the thick torus. Even if that is the case, our result is 
consistent with the abundance ratio. 

In this regime of e = 0.3 - 0.9 and M BH = 1O 7_8 - 5 M , 
our model indicates that the observer's viewing angle is 
responsible for the difference between BAL and non-BAL 
quasars. If the observer is in the direction of the wind 
flow the absorption features are detected. Otherwise, the 
nuclei arc identified as non-BAL quasars. That is, our re- 
sult gives physical bases in phenomcnological model (Elvis 
2000), in which the dichotomy between the BAL and non- 
BAL quasars is explained by the observer's viewing angle. 
Note that, in addition to such a viewing-angle effect, our 
results imply that the BAL features does not appear if 
the Eddington ratio is too small (e < 0.01) or if the black 
holes is too massive (Mbh ^ 10 9 M Q ). 

Here, we note that our model overestimate the BAL 
probability in the regime of e = 0.05 — 0.1. Although our 
model gives the probability of 20 — 50%, the accurate value 
would be much small. The launching matter is blown away 
for a wide angle for e = 0.05 — 0.1, leading to the larger 
BAL probability. The streamlines for e = 0.1 is plotted 
in Figure 6. In this case, the BAL probability is 48 %. 
However, then, the mass outflow rate of the wind is a 
few times larger than the mass accretion rate of the disk. 
This implies that our model overproduces the outflows, 
and actual BAL probability should be smaller than 48 %. 
In the case that the mass accretion rate is comparable to 
the mass outflow rate, we need to solve the outflow and 
the disk, self-consistently (we will discuss later). 

4-2. Comparison with previous works 

Throughout the present study, we employ the non- 
hydrodynamic method that is basically the same with that 
used in Risaliti & Elvis (2010). In Risaliti & Elvis (2010), 
the initial density and velocity are treated as free param- 
eters and assumed to be constant in the radial direction, 
for simplicity. However, we employ more realistic profiles 
of po and vq based on CAK75, in which the mass outflow 
rate of the line-driven stellar wind is investigated in detail. 
Such a discrepancy of the initial values leads to the dif- 
ference of the wind structure. In Risaliti & Elvis (2010), 
since the matter is overionized by X-ray irradiation at the 
vicinity of the black hole, the line force is made power- 
less and fails to launch the wind. Murray et al. (1995) 
also reported that the failed wind appears in the vicinity 
of the black hole via X-ray irradiation. In our work, the 
matter is not highly ionized via the high density (see po 
profiles in Figure 2). However, high density reduces the 
force multiplier and prevents the launching of the wind as 
we have mentioned in section §3.1. Also, the outer part 
of the disk contributes more to the total mass loss rate 
(oc poVar 2 ) than the inner disk in Risaliti & Elvis (2010), 
since po and vq are constant in the radial direction. In 
contrast, po as well as Vq increase with the decrease of r, 
meaning that the disk loses more mass from the inner part 
of the launching region. 
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Fig. 6. Same as Figure 3, but for e = 0.1. In this case, the 
wind with large covering factor forms. 

Here, we note that the wind structure calculated by 
the present non-hydrodynamic method is similar to that 
calculated by hydrodynamic simulations of PK04. The 
opening angle of the wind, ~ 60°, for PK04 approximately 
coincides with our result, ~ 50°. The column density of 
the wind in the moderate or low ionization state is around 

10 24-25 cm -2 (fl _ 70 o) in p K Q 4j and 10 24-25 cm -2 (0 _ 

60°) and > 10 25 cm" 2 (0 > 70°) in our model. The wind 
velocity is much larger in our work than in PK04. The 
outward velocity is 100,000kms _ in our work, but PK04, 
on the other hand, showed around 10,000kms _1 (the peak 
value is 20, 000 km s" 1 ). Such a discrepancy is though to 
be caused by that the disk wind is launched from more 
inner region in our model than in PK04. In PK04, as 
well as Risaliti & Elvis (2010), the density just above the 
disk surface is set to be constant in the radial direction. 
Thus, the density at the wind base is lower in PK04 than 
in our model at the inner region. The less dense matter is 
almost fully ionized and is not accelerated by line force. 
The failed wind region in PK04 is larger than that in 
our model. The wind launched from the inner region is 
effectively accelerated by strong UV radiation from the 
inner part of disk, leading to higher velocity of the wind. 
Here we comment that, when the initial density is set 
to be the same as that in PK04, we can approximately 
reproduce the results of PK04, e.g., the wind velocity of 
10,000kms _1 and the failed wind due to overionization. 

Feldmeier et al. (1999) investigated the line driven wind 
in cataclysmic variables. Their streamlines are almost in 
pararell with each other. Such a feature is similar with out 
resutls. However, they reported that the opening angle of 
the wind is 10° — 30°, which is much smaller than that 
of our model and PK04, 50° — 60°. Although the stream- 
lines bends above the disk in our model, straight lines are 
assumed in Feldmeier et al. (1999). Also, the dilution of 
the radiative flux by the wind (self-shielding) is not taken 
into consideration in Feldmeier et al. (1999). 



4-3. Future work 

Although axisymmctric hydrodynamic simulations of 
the line-driven wind are attempted by PK00 and PK04, 
three-dimensional hydrodynamic approaches would be im- 
portant future work. Recently, the time variation of the 
absorption features has been reported by Misawa et al. 
(2007). This implies that the wind structure changes with 
time and/or non-axisymmetric pattern exists in the winds. 
Such problems would be resolved by three-dimensional 
hydrodynamic simulations. Owocki & Puis (1999) have 
reported using one-dimensional calculations that the line- 
driven wind exhibits violent density fluctuations and many 
density peaks. In addition, Proga et al. (2000) and 
PK04 have shown that high density blobs are sometimes 
generated in the wind. However, their work is a two- 
dimensional study. We will explore the three-dimensional 
simulations to investigate the time evolution and non- 
axisymmetric structure of the winds. 

In the case that the mass outflow rate is comparable 
to the mass accretion rate of the disk, we should solve 
the wind and disk, self-consistently. The wind extract the 
mass, angular momentum, and the energy from the disk. 
The disk structure, thus, changes via launching of the 
wind, and the wind structure is subject to influence of the 
changing the disk structure. As we have discussed above, 
the mass outflow rate is comparable to or larger than mass 
accretion rate for e = 0.05 — 0.1 in our model. By X- 
ray observations of AGNs, Tombesi et al. (2012) reported 
that the mass outflow rate via the ultra fast outflows is 
comparable to the mass accretion rate onto the black hole. 

In the present study, we verify the line-driven wind 
model based on the ionization parameter, the outward ve- 
locity, and the column density. However, we should calcu- 
late the emergent spectra and directly compare these with 
the observations. Such a study was recently attempted by 
Sim et al. (2010), in which they performed Monte Carlo 
radiative transfer simulations in the X-ray band. 

5. Conclusions 

We have studied the structure of the disk wind driven 
by the radiation force including the line force, by calculat- 
ing trajectories of the fluid elements which are launched 
from the surface of the geometrically thin and optically 
thick disk. Here, the density and velocity at the wind base 
are set so as meet the condition of line-drive stellar wind 
(CAK75). We have solved the equation of motion, cou- 
pling with the mass conservation along the streamlines. 
The radiation force, the ionization parameter, and the 
temperature are calculated by taking into consideration 
the extinction of the X-ray from the vicinity of the black 
hole and the UV from the disk by the ejected matter. We 
have compared the resulting wind structure with the wind 
properties of X-ray observations of the BAL quasars (the 
ionization parameter, the outward velocity, and the col- 
umn density), and estimated the probability (BAL prob- 
ability), with which the system is identified as a BAL 
quasar. 
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In the case that the Eddington ratio is e ~ 0.3 — 0.9 
and the black hole mass is Mbh ~ 1O 7_8 ' 5 M0, we found 
that the funnel-shaped disk wind with a opening angle of 
~ 50° forms. In this regime, since the wind shape is insen- 
sitive to the Eddington ratio and the black hole mass, the 
BAL probability of several to 10 percent does not change 
so much. The observer's viewing angle is responsible for 
whether or not the quasars exhibit the broad absorption 
lines. Our model is consistent with the phenomenological 
model of BAL quasars proposed by Elvis (2000). 

The wind is launched from the middle region of the 
disk, where the line force is strong enough to accelerate 
the matter. At the inner region, since the larger initial 
density reduces the force multiplier (line force), the ejected 
matter immediately returns to the disk surface. On the 
other hand, in the outer region, since the UV radiation 
is attenuated via the obscuration by the matter in the 
inner and middle regions, the radiation force fail to launch 
the wind. In the case of Mbh = lO 7 Af and e = 0.3, the 
launching region is 30i?s ^ r 350i?s- 

If the Eddington ratio is very small (e < 0.01), the disk 
wind does not appear, since the UV intensity is too week 
for the matter to be accelerated by the line force. Also, 
the disk wind is not launched when the black hole is too 
massive (Mbh ^ 10 9 A/©). This is because that the force 
multiplier at the wind base decreases with an increase of 
the black hole mass. In addition, the disk temperature 
also decreases (UV emission region, Toff > 10 4 K, shrinks) 
with an increase of the black hole mass. Thus, for e < 0.01 
or for Mbh ^ 1O 9 M , the BAL quasars are not observed 
independently of the viewing angle. 
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